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Executive  Summary 


The  Air  Force  has  requirements  for  accurate  knowledge  and  specification  of  the 
state  variables  of  the  Earth's  thermosphere.  Work  done  under  this  contract  will  lead  to  a 
new  model  of  global  thermospheric  density  that  can  be  used  to  specify  and  forecast  neutral 
densities,  temperatures,  and  winds  for  a  wide  range  of  solar  and  geomagnetic  activity.  The 
model  will  be  based  on  simulations  made  with  the  National  Center  for  Atmospheric 
Research  (NCAR)  Thermosphere/Ionosphere  General  Circulation  Model  (TIGCM)  and  on 
data.  It  will  be  capable  of  using  real-time  geophysical  indices  or  data  from  grmmd-based 
and  satellite  inputs  and  will  provide  neutral  variables  at  specified  locations  and  times  in  the 
altitude  range  90  -  1500  km.  This  hybrid  theoretical  and  semi-empirical  model  will  be  based 
on  a  new'  Vector  Spherical  Harmonic  (VSH)  analysis  technique  developed  at  the  University 
of  Michigan  that  permits  the  incorporation  of  the  TIGCM  outputs  and  data  into  the  model. 
The  VSH  model  will  be  provided  for  use  by  the  Air  Force  in  an  operational  setting  and  will 
be  a  more  accurate  version  of  existing  models  of  the  neutral  upper  atmosphere. 

This  technical  report  describes  the  progress  of  the  University  of  Michigan  and 
National  Center  for  Atmospheric  Research  teams  in  the  development  of  the  hybrid  model. 
Scientific  progress  and  plans  are  reported  for  each  quarterly  reporting  interval  throughout 
the  first  year  of  the  project. 
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The  Air  Force  has  requirements  for  accurate  knowledge  and  specification  of  the 
state  variables  of  the  Earth's  thermosphere.  Work  done  under  this  contract  will  iead  to  a 
new-  model  of  global  thermospheric  density  that  can  be  used  to  specify  and  forecast  neutral 
densities,  temperatures,  and  winds  for  a  wide  range  of  solar  and  geomagnetic  activity.  The 
model  will  be  based  on  simulations  made  with  the  National  Centc  for  Atmospheric 
Research  (NCAR)  Thermospheric/Tonospheric  General  Circulation  Model  (TIGCM)  and  on 
data.  The  modei  will  be  capable  of  using  real-time  geophysical  indices  or  data  from  ground- 
based  and  satellite  inputs  and  will  provide  neutral  variables  at  specified  locations  and  times 
in  the  altitude  range  90  -  1500  km.  The  model  will  be  semi-empirical  and  will  be  based  on  a 
new  vector  spherical  harmonic  (VSH)  analysis  technique  developed  at  the  University  of 
Michigan  that  permits  the  incorporation  of  the  TIGCM  outputs  into  the  model.  The  model 
will  be  provided  for  use  by  the  Air  Force  in  an  operational  setting  as  a  more  accurate 
version  of  existing  models  of  the  neutral  upper  atmosphere. 


Plans  for  the  next  quarter  of  effort  relate  to  the  work  statement  given  in  the  original 
proposal.  The  timeline  for  model  development  is  given  in  Appendix  1.1.  Work  to  verify  the 
TIGCM  in  the  altitude  region  above  250  km  is  underway.  Two  major  case  studies  are 
planned,  one  involving  10  days  of  data  from  the  Dynamics  Explorer  satellite  in  198?  and 
one  involving  a  previous  isolated  geomagnetic  storm  study.  The  Dynamics  Explorer  data 
has  been  provided  to  Dr.  Roble  and  a  series  of  TIGCM  runs  are  planned  that  will  be 
extensively  analyzed  as  part  of  the  verification  effort.  Dr.  Roble  has  prepared  a  computer 
animation  of  the  isolated  storm  study  that  will  be  shown  at  the  CEDAR  meeting  in  June 
1989.  Considerable  effort  at  the  University  of  Michigan  will  be  made  to  speed  up  the 
current  VSH  algorithm.  It  is  anticipated  that  an  increase  in  speed  of  about  a  factor  of  6-10 
can  be  realized.  Further  work  will  be  carried  out  to  define  the  VSH  library  approach.  Work 
will  be  also  carried  out  to  define  the  TIGCM  input  stategy  for  this  program,  as  well  as  to 
incorporate  data  from  Dynamics  Explorer  into  the  VSH  model  for  a  few  cases  using  an 
objective  analysis  scheme. 
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Executive  Summary. 


The  Air  Force  has  requirements  for  accurate  knowledge  and  specification  of  the 
state  variables  of  the  Earth's  thermosphere.  Work  done  under  this  contract  will  lead  to  a 
new  model  of  global  thermospheric  density  that  can  be  used  to  specify  and  forecast  neutral 
densities,  temperatures,  and  winds  for  a  wide  range  of  solar  and  geomagnetic  activity.  The 
model  will  be  based  on  simulations  made  with  the  National  Center  for  Atmospheric 
Research  (NCAR)  Thermospheric/Ionospheric  General  Circulation  Model  (TTGCM)  and  on 
data.  The  model  will  be  capable  of  using  real-time  geophysical  indices  or  data  from  ground- 
based  and  satellite  inputs  and  will  provide  neutral  variables  at  specified  locations  and  times 
in  the  altitude  range  90  -  1500  km.  The  model  will  be  semi-empirical  and  will  be  based  on  a 
new  vector  spherical  harmonic  (VSH)  analysis  technique  developed  at  the  University  of 
Michigan  that  permits  the  incorporation  of  the  TIGCM  outputs  into  the  model.  The  model 
will  be  provided  for  use  by  the  Air  Force  in  an  operational  setting  as  a  more  accurate 
version  of  existing  models  of  the  neutral  upper  atmosphere. 

During  the  past  reporting  period,  the  effort  at  the  University  of  Michgian  has 
focussed  on  the  development  and  improvement  of  the  basic  algorithm,  based  on  the  vector 
spherical  harmonic  expansion  of  TIGCM  calculations.  The  speed  of  the  routine  has  been 
increased  by  a  factor  of  6  by  developing  more  efficient  loops,  etc.  In  addition, 
documantation  has  been  produced  to  facilitate  testing  of  the  routine  by  AFGL.  A 
preliminary  version  of  the  VSH  model  was  sent  to  AFGL  for  evaluation,  complete  with  12 
sets  of  coefficients,  documentation  and  source  code.  In  addition,  a  perliminary  version  of 
the  planned  objective  analysis  scheme  for  the  ingestion  of  experimental  data  was 
developed,  using  thermospheric  wind  measurements  from  Dynamics  Explorer  as  a  first  test 
case.  This  work  led  successfully  to  the  development  of  a  semi-empirical  model  of  polar 
thermosphere  dynamics. 

Work  at  the  National  Center  for  Atmospheric  Research  has  continued  to  focus  on 
improvements  to  the  TIGCM  and  on  validation  studies  involving  data/theory  comparisons. 
A  specific  comparison  between  TIGCM  calculations  and  SETA  data  has  produced 
encouraging  results  indicating  that  the  TIGCM  can  predict  density  perturbations  with 
reasonable  accuracies  (see  below  for  details). 
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Scientists  from  the  U nviersity  of  Michigan  and  from  NCAR  attended  the  second 
quarterly  review  meeting  in  Boulder  on  21  September  1989.  This  meeting  was  also 
attended  by  many  other  scientists  from  various  institutions  as  well  as  Air  Force  Personnel. 

Plans  for  the  next  reporting  period  include  the  further  development  of  the  VSH 
approach,  leading  to  delivery  of  the  Phase  1  Model  at  the  end  of  FY  1990.  Additional  key 
validation  and  model  improvement  efforts  will  be  carried  out  at  NCAR. 

2.1  Progress  during  therurrent  reporting  period. 

Progress  during  the  current  reporting  period  was  presented  at  the  second  quarterly 
meeting  in  Boulder,  Colorado  on  September  21  1989.  The  detailed  agenda  for  this  meeting 
is  shown  in  Appendix  2.1.  Presentations  were  made  by  Prof  Killeen  of  the  Unviersity  of 
Michigan  and  Dr.  Roble  of  NCAR  and  materials  from  their  presentations  are  included  as 
appendicies  to  this  report. 

Prof.  Killeen  presented  the  work  carried  out  in  association  with  the  VSH 
development  effort.  Appendix  2.2  presents  current  documentation  for  the  model,  including 
a  "top  down"  structural  design  description,  a  list  of  available  coefficient  sets,  a  model 
variable  directory,  and  a  description  of  the  input  and  output  variables. 

Prof.  Killeen  also  discussed  the  schedule  for  the  effort  (see  Appendix  2.3).  The 
overall  model  development  effort  is  broken  down  into  13  sub-tasks.  Tasks  1-5  have  been 
initiated  and  task  1  completed.  Although  the  program  had  a  late  start  in  March  1989,  the 
work  is  on  schedule  for  delivery  of  the  Pase  I  model  at  the  end  of  FY  1990. 

The  University  of  Michigan  subcontract  to  NCAR  has  been  established  and  a 
division  of  responsibilities  for  the  team  at  the  two  institutions  has  been  agreed  upon.  An 
initial  bench-test  version  of  the  VSH  model  has  been  delivered  to  AFGL  for  evaluation.  It 
is  anticipated  that  early  iterations  of  this  type  will  ensure  an  improved  Phase  1  Model. 

An  objective  analysis  scheme  was  developed  for  the  incorporation  of  neutral  wind 
measurements  from  Dynamics  Explorer  into  the  VSH  model.  Averaged  wind 
measurements  for  the  high  latitude  regions  of  both  hemispheres  were  binned  according  to 
IMF  By  and  Kp  indices.  The  averaged  winds  were  merged  together  with  the  TIGCM 
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predictions  for  corresponding  geophysical  conditions  and  converted  to  a  form  for  inclusion 
in  the  VSH  model.  The  coefficients  sets  for  these  cases  were  delivered  to  AFGL  for 
evaluation. 

A  total  of  106  TIGCM  runs  have  been  identified  as  necessary  for  the  Phase  1  Model 
development  (see  Appendix  2.3).  These  runs  will  be  made  during  the  currei.i  fiscal  year. 

Dr.  Roble  reported  on  the  progress  made  in  updating  and  performing  validation 
studies  on  the  TIGCM  (see  Appendix  2.4).  One  of  the  key  case  studies  has  involved  the 
detailed  one-on-one  comparison  with  total  density  data  along  the  track  of  the  Air  Force 
satellite  containing  the  SETA  instrument.  The  model  was  tuned  specifically  to  match  the 
time  dependent  thermospheric  forcings  corresponding  to  the  March  22-23  1979  data  set. 
Appendix  2.4  shows  the  results  of  this  study.  Percent  differences  for  measured  minus 
model  densities  are  shown,  indicating  a  standard  deviation  about  the  mean  of  ~9%  for 
day  side  data  and  -12%  for  night  side  data.  It  should  be  noted  that  the  TIGCM  total 
densities  were  corrected  by  a  constant  multipicative  factor  in  order  to  match  the  measured 
DC  level.  This  correction  factor  implies  that  the  total  heating  rate  in  the  model  still  needs 
refinement  to  match  the  observed  total  density  structure.  Thus  work  is  ongoing  to  tune  the 

various  aeronomical  heating  and  cooling  rates  to  match  better  the  measurements. 

■> 

Nevertheless,  the  agreement  with  between  the  modelled  and  measured  density  variations 
about  the  global  mean  are  extremely  encouraging  and  indicate  that  the  TIGCM  is  already 
providing  an  improvement  over  pure  semi-empirical  models,  such  as  MSIS. 
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22  Flans  for  the  next  reporting  period 


Plans  for  the  next  quarter  of  effort  relate  to  the  work  statement  given  in  the  original 
proposal,  see  Appendix  2.3.  Specifically,  work  on  tasks  2-5  will  continue. 

Work  to  verify  the  TIGCM  in  the  altitude  region  above  250  km  will  continue. 
Specific  tests  involve  detailed  comparisons  with  satellite  data  from  a  ten  day  period  (tidy.-. 
83021  -  83031)  near  the  end  of  the  Dynamics  Explorer  mission  when  full  coverage  of  both 
hemispheres  was  possible  with  all  instruments  functioning. 

Two  major  case  studies  are  planned,  one  involving  the  10  days  of  data  from  the 
Dynamics  Explorer  satellite  in  1983  and  one  involving  a  previous  isolated  geomagnetic 
storm  study.  The  Dynamics  Explorer  data  has  been  provided  to  Dr.  Roble  and  a  series  of 
TIGCM  runs  are  planned  that  will  be  extensively  analyzed  as  part  of  the  verification  effort. 
Dr.  Roble  has  prepared  a  computer  animation  of  the  isolated  storm  study  that  will  be  shown 
at  the  next  quarterly  meeting  in  January.  Further  work  will  be  carried  out  to  define  the 
VSH  library  approach.  Work  will  be  also  carried  out  to  define  the  TIGCM  input  stategy  for 
this  program,  as  well  as  to  incorporate  data  from  Dynamics  Explorer  into  the  VSH  model 
for  a  few  cases  using  an  objective  analysis  scheme. 
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VSH  TOP-DOWN  DESIGN 


Upon  tich  call  to  tha  VSH  pcograa,  tha  following  eight  atap  procedure  is 
carried  out.  Tha  pciaary  purpose  of  aach  step  and  tha  subcoutina  responsible 
for  executing  tha  step  is  given  below. 

Each  of  tha  indicated  subroutines  is  called  froa  tha  aain  controlling 
routine  VALUE.  A  aain  driving  program  is  then  needed  to  call  VALUE. 

1.  Detaraine  whether  any  of  the  following  have  changed  since  last  call: 

Latitude,  tiae,  TIGCM  run,  variables  requested  I  INITIAL) 

2.  Detaraine  if  all  needed  variables  have  been  requested? 

Densities  need  teaperatura  and  aixing  ratios. 

Zonal  and  meridional  winds  need  each  other.  (  INITIAL) 

J.  If  TIGCM  run  or  variables  requested  have  changed: 

Read  coefficients  for  requested  variables  of  desired  TIGCM  run 
site  is  (9  field  variables,  3  altitudes,  5  tonal  wavenumbers, 

IS  meridional  wavenuabar a,  7  Fourier  tiae  coefficients). 

4.  If  tiae,  run,  or  variables  have  changed: 

Evaluate  Fourier  series  for  the  requested  tiae  at  each  variable,  altitude, 
tonal  wavenuabar,  and  aeridional  wavenuabar.  Coefficient  array  site  is 
now  (9,3,5,251.  (TIHETR) 

5.  If  latitude  has  changed : 

For  the  requested  latitude,  find  the  values  of  the  vector  spherical 
haraonics  P,  v,  and  W  for  each  tonal  and  aeridional  wavenuabar.  IBAS1S) 

6.  Evaluate  the  VSH  representation  at  the  requested  latitude  and  longitude  for 
each  variable  at  each  altitude: 

For  neutral  and  ion  winds,  usa  vector  spherical  haraonics  V  end  w. 

For  other  variables,  use  scalar  spherical  harmonic  P. 

Coefficient  array  site  la  now  (9,3).  (SPACTR) 

7.  At  each  altitude,  derive  densities  O,  N2,  and  02. 

Coefficient  array  site  is  now  (12.3)  (DENSITY) 

8.  If  a  specific  altitude  was  raquestad  (rather  than  a  pressura  surfaca), 
in t a r po 1  a t ion  or  extrapolation  is  necessary: 

For  taaperature,  asauae  a  Bates  fit. 

For  densities,  essuae  a  linear  fit  in  their  logarithms. 

For  other  variables,  assuae  a  quadratic  fit. 

However,  if  e x t r a po  1  a t i on  is  needed  above  the  highest  pressure  surface, 
assuae  that  winds  and  taaparature  have  reached  their  asymptotic  values. 

I ALTTR ) 


MODIFICATIONS  IN  PROGRESS 

The  vertical  interpolation  procedure  (step  8)  will  be  changed  into  a 
Chebyshev  function  representation.  This  eodification  will  allow  a  more 
accurata  representation  in  tha  vertical  of  tidal  affects. 

The  Fourier  representation  in  time  (step  J  )  will  also  be  changed  ir. to  a 
Chebyshev  basis.  This  modification  will  allow  a  more  accurate  representation 
of  magnetic  stora  conditions  las  the  restriction  to  periodic  boundary  conditions 
will  be  eliminated  . 


Array 
(  VSHIN  ) 
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TtGCM  RUNS  AVAILABLE  TOR  VSH 


1  . 


2  . 


3  . 


4  . 


5  . 


6  . 


7  . 


8  . 


9  . 


1  0  . 


1  1  . 


I  2  . 

1  3  . 

I  }  . 


Solar  aaxiaua  winter  solstice  conditions  with  average  geomagnetic 
activity  (Kp=3.  Ap  20,  F10. 9=220,  Croas-cap  potential  60kV) 

Solar  ainiaua  winter  solstice  conditions  with  average  geoaagnetic 
activity  (Kp=3.  Ap  20,  F10. 7=120,  Cross-cap  potential  30kV) 

♦ 

Solar  liniaua  aquinoi  conditions  with  average  geoaagnetic  activity 
(Kp=3,  Ap  20,  F10. 7=120,  Cross-cap  potential  30kV) 

Solar  aaxiaua  solstice  conditions  with  average  aagnetic  activity 
By  positive  conditions 

Solar  aaxiaua  solstice  conditions  with  average  aagnetic  activity 
By  negative  conditions 

Model  transforaed  into  aagnetic  coords,  for  solar  aaxiaua  at  solstice 
with  average  aagnetic  activity  By  negative  conditions 

Model  transforaed  into  aagnetic  coords,  for  solar  aaxiaua  at  solstice 
with  average  aagnetic  activity  By  positive  conditions 

Solar  aaxiaua,  Dec  -solstice  avg  activity  with  B Y -0 , hp» 1 1  , x -c a p  =  6 0 k V  is 
aerqed  with  DE-2  data.  The  DE-2  winds  are  for  Kp  < a  3.0,  By  >  0.0  with 
the  h l g h -  1  a t l t ude  data  (  >  3  5  degrees)  in  geoaagnetic  coords,  and 
the  low  latitude  data  in  geographic. 

Solar  aaxiaua,  Dec  solstice  avg  activity  with  BY  =  0  ,  hp= 1 1 , x -cap  =  6 0 k V  is 
aerged  with  DE-2  data.  The  DE-2  winds  are  for  Kp  >  3.0,  By  >  0.0  with 

the  h i gh -  1  a 1 1 t ude  data  (>35  degrees)  in  geoaagnetic  coords,  and 

the  low  latitude  data  in  geographic. 

Solar  aaxiaua,  Dec  solstice  avg  activity  with  BY=»0  ,  hp«  1 1  ,  x -cap=>6  0k  V  is 
aerged  with  DE-2  data.  The  DE-2  winds  are  for  Kp  <=  3.0,  By  <  0.0  with 

the  h i gh- 1  a t i t ude  data  <>35  degrees)  in  geoaagnetic  coords,  and 
the  low  latitude  data  in  geographic. 

Solar  maximum,  Dec  solstice  avg  activity  with  B  Y  =>0  ,  hp=  1 1  ,  x -ca  p  =  6  0  k  V  is 

aerged  with  DE-2  data.  The  DE-2  winds  are  for  Kp  >  3-0,  By  <  0.0  with 

the  high-latitude  data  (>35  degrees)  in  geomagnetic  coords,  and 

the  low  latitude  data  in  geographic. 

Solar  aaxiaua  solstice  conditions  with  low  geoaagnetic  activity 
( K  p  = l  ,  H  P  =  3  ,  F10.  7  =  220,  Cross-cap  potential  30kV) 

Solar  a  a  x  l  r,  ui  solstice  condition  •••ith  average  geomagnetic  activity 
(Kp=3,  H  P  ■  11,  F10. 7-220,  Cross-cap  potential  60kV) 

Solar  maximum  solstice  conditions  with  average  geomagnetic  activity 
(  K  p  =  5  ,  H  P  =  33,  F10. 7=220,  Cross-cap  potential  90kV( 
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VSH  VARIABLE  DICTIONARY 


ALT  REAL  Requested  altitude  or  pressure  surface  (t,  2,  or  3) 

FLDDIM  INTEGER  Dimensioned  number  of  variables  (excluding  densities) 

FLDMAX  INTEGER  Number  of  variables  actually  used  (excluding  densities) 

GCOLAT  REAL  Requested  geographical  colatitude 

GLON  REAL  Requested  geographical  longitude 

GETDENS  LOGICAL  -TRUE,  if  a  call  to  DENSITY  is  needed 

HRDIM  INTEGER  Dimensioned  number  of  time  periods  in  coefficient  data 

HRMAX  INTEGER  Number  of  time  periods  represented  in  coefficient  data 

I  PR  I NT  INTEGER  A  debug  print  control;  not  presently  in  use 

IRUN  INTEGER  Requested  TIGCH  run  number 

LASTLAT  REAL  Requested  latitude  on  previous  call 

LASTRUN  INTEGER  Requested  TIGCH  run  on  previous  call 

LASTUSI  LOG I  CAL  (  12  )  Requested  USING  values  on  previous  call 

LASTTIN  REAL  Requested  UT  on  previous  call 

MDIM  INTEGER  Dimensioned  number  of  zonal  wavenumbers 

UMAX  INTEGER  Number  of  zonal  wavenumbers  actually  used  this  run 

NEWLAT  LOGICAL  .TRUE,  if  a  new  latitude  has  been  requested 

NEWTIM  LOGICAL  .TRUE,  if  a  new  UT  has  been  requested 

NEWRUN  LOGICAL  .TRUE,  if  a  new  run  has  been  requested 

ONES  INTEGER  IRUN  modulo  10 

PBAR  RE  AL  (  HDIH-t  1  ,  NDI M+- 1  )  Scalar  spherical  harmonic  P 

PHASE  REAL  Phase  angle 

SOLAR  INTEGER  0  or  1  depending  on  solar  max  or  solar  min 

STVSH  REAL! ZDIN,  VARDIH  I  Coefficient  array  following  space  and  time 

transformation 

TDIH  INTEGER  Dimensioned  number  of  time  periods 

TENS  INTEGER  I n t e g e r ( I RUN/ 1 0  I 

T VS  H  REAL)  FLDOIH,  ZD  IH  , MDIM, NDIM)  Coefficient  array  following  time 

transformation 

TMAX  INTEGER  Number  of  time  periods  actually  used 

USING  LOGICAL! 12  )  .TRUE,  for  a  variable  if  it  has  been  requested 

UT  REAL  Requested  universal  time 

UTHR  REAL(HRDIM)  UT  '  s  specified  in  coefficeint  file 

VARDIM  INTEGER  Dimensioned  number  of  variables  including  densities 

VARMAX  INTEGER  Number  of  variables  including  densities  actually  used 

VBAR  REAL! MDIM+l  ,  NDIM  +  l  )  Vector  spherical  harmonic  V 

VSH  REAL  (  FLDDIM . ZDIM , MDIM , NDIM , TDIH )  Coefficient  array  as  read  in 

from  data  files 

XALT  REAL(ZDIM)  Altitudes  specified  in  coefficient  file 

WBAR  REAL  (  MD IM+ 1 , ND I M+ 1 )  Vector  spherical  harmonic  W 
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COEFFICIENT  DATA  FILES 


Currently,  14  runs  of  the  TIGCM  model  are  available  to  the  VSH  program, 
ach  is  stored  in  an  identical  binary  unformatted  form  in  files  RUN01.DAT 
Trough  RUN14.DAT.  For  each  of  the  9  field  variables,  there  are  two  records: 
head°r  record  and  the  data  record. 

The  header  record  contains  information  on  the  truncation  level  of  the 
oefficient  data  as  well  as  their  altitudes  and  universal  times.  The  data  is 
tored  in  the  sequence:  MM AX ,  NMAX,  TMAX ,  ZMAX,  HRMAX ,  (UT11R(I),  I=1,HRMAX), 
XALT(I),  1=1, ZMAX). 

The  data  record  is  stored  as  ( ( ( ( VSH ( VAR , Z ,M, N , T) ,  Z=1,ZMAX),  M-l.MMAX), 
=  1 , NMAX) .  M= 1 , MMAX  > 
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VSR  INPUT/OUTPUT  VARIABLES 


Nik 

OUTPUT ( 1 2 ) 


Ni«l 

USING! 12  ) 

GLON 

GCOLAT 

UT 

ALT 

IRUN 


OUTPUT  VARIABLES 


Type  Description 


REAL  the  array  of  output  variables: 

tonal  neutral  wind  (U) 

Meridional  neutral  wind  (V) 
neutral  temperature  (T) 
meridional  ion  drift  f(Ui) 
tonal  ion  drift  ( Vi ) 
vertical  neutral  wind  (W) 

height  of  the  content  pressure  level  (ALT) 
atomic  oxygen  mass  mixing  ratio 
Molecular  nitrogen  mass  nixing  ratio 
atomic  oxygen  number  density 
molecular  oxygen  number  density 
molecular  nitrogen  number  density 


INPUT  VARIABLES 


Typo  Description 


LOGICAL 


REAL 

REAL 

REAL 

REAL 


INTEGER 


an  array  of  12  logical  variables  (.true,  or  .false.) 
corresponding  to  each  of  the  above  12  variables; 
setting  a  value  to  .false,  eliminates  calculation  of 
of  that  variable  (to  reduce  computational  time) 
geographic  longitude  GLON  (degrees) 
geographic  co-latitude  GCOLAT  (degrees) 

Universal  Time  UT  (hours) 

altitude  (km)  (between  130  and  600  km);  or 

a  pressure  surface  (indicated  by  a  1,  2,  or  3): 

ALT*  l:  denotes  -4  surface  ("120km) 

ALT*  2:  denotes  z*  -l  "  ('185km) 

ALT*  3:  denotes  r*  +1  "  |'350km) 

(For  solar  Minimum  cases,  z*-4,  -2,  and  +2) 

the  particular  TIGCM  run: 
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NTM  -  AFGUAWS,  Boulder  5/8/4  TLK  1 


26 


NTM  -  AFGUAWS,  Boulder  5/8/4  TLK  3 
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program.  Status  of  VSH/Objective  analysis  algorithms,  TIGCM 
Validation,  etc. 


NTM  -  AFGUAWS,  Boulder  5/8/4  TLK  9 


Phase  I  TIGCM  runs 
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NTM  -  AFGL/AWS,  Boulder  5/8/4  TLK  10 
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NCAR  THERMOSPHERE/IONOSPHERE 
GENERAL  CIRCULATION  MODEL  (TIGCM) 


Primitive  equations  of  dynamic  meteorology  adapted  to 
thermospheric  heights. 

Eulerian  model  of  the  ionosphere. 

Horizontal  grid  5  *  latitude  X  5  °  longitude,  geo¬ 
graphic. 

Vertical  grid  -  25  constant  pressure  surfaces,  2  grid 
points  per  scale  height,  95  to  500  km. 

Time  step  150  or  300  S. 


INPUT 


Solar  EUV  and  UV  radiation  5  to  250  nm. 

Empirical  ionospheric  convection  and  auroral  particle 
precipitation  models. 

Structure  of  upward  propagating  tides  from  middle 
atmosphere. 
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Fig.  6 
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Fig.  7 
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sec-1) 


TGCM  NEUTRAL  TEMPERATURE  (K)  *  (U.V)  (M/S) 
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YD  =  76080  AVE  HT  =  279.5 
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LATITUDE  ,  .  LATITUDE 
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HEELIS  VI  (M/5)  ♦  POTENTIAL  (KEV) 

LOCAL  TIME 

ECR245.246  12 

YD  =  76080 

UT  =  19.00  /  9- 
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BY  = 
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J  f 


FRAME  301 


800  M/S 


CONTOUR  TROW  -20.000  TO  20.000  CONTOUR  INTERVAL  Of  >0.000  RT(J,3)«  1.J047 


Trriy(  Kjn  ijpai  TEMPERATURE  (K)  ♦  (U.V)  (M/S)  27.5  N 

LOCAL  TIME 

EQSN03.N04  12  Z  -  2. 
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TAIL  POTENTIAL  KV 


1DIAL  NO.  POINlb*  7766 


ALL  ff4€M  OAT  SKJE  DATA  TOTAL  1C.  POINTS-  7766 
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HERCEHT  DIFFERENCE  FROM  MSIS  B6 

Fig.  20 


50 


OF  TOTAL  POINTS  PERCENT  OF  TOTAL  POINTS 


ALL  N-FEM  NIGHT  SICE  DATA 

1 2  i - 1 - 1 - 1 - 1 - 1 - 1  i 


TOTAL  NO.  POINT S=  8140 


PERCENT  DIFFERENCE  FROM  TIGCM  (NORMALIZED  BY  0.75B0) 


ALL  N-FEM  NIGHT  SIDE  DATA  TOTAL  NO.  POINT S=  8140 


PERCENT  DIFFERENCE  FROM  MSIS  06 

Fig.  ?1 
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Croes-cAp  Potential  (calc  fro* 
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ut  tmi 


OF  TOTAL  POINTS  PERCENT  OF  TOTAL  POINTS 


PERCENT  DIFFERENCE  FROM  ftSIS  86 

Fin.  25 
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00 
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SUMMARY  AND  CONCLUSIONS 


•  TIGCM  used  to  simulate  time-dependent  variations 
frocn  20  March  to  29  March  1979. 

•  Auroral  particle  precipitation  specified  by  Hemisphere 
Power  Index. 

•  Ionospheric  convection  pattern  variations  prescribed  by 
solar  wind  variations. 

•  Results  show  considerable  variability  in  all  TIGCM 
fields  in  response  to  high  latitude  forcings,  e.g.,  latitu¬ 
dinal  surges  and  recovery. 

•  Test  of  aeronomic  scheme  stability  and  ability  to  relax 
to  pre-storm  conditions. 

•  Need  to  develop  post  processors  for  analysis  of  physi¬ 
cal,  chemical  and  radiative  processes,  e.g.,  wave  energy, 
energy  of  mean  flow,  energy  radiated  etc. 
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SUMMARY  AND  CONCLUSIONS 


•  TIGCM  compares  reasonably  well  with  overall  varia¬ 
tions  predicted  by  MSIS-86,  but  with  considerably 
more  structure  and  variability. 


•  Compared  SETA-1  measured  density  with  MSIS-86  and 
TIGCM  predictions  along  satellite  path  for  the  9  day 
period: 

•  TIGCM  and  MSIS-86  give  comparable  agreement 
with  dayside  measurements. 

•  MSIS-86  performs  better  at  night  than  TIGCM. 


•  Need  to  improve  plasma  exchange  processes  between 
magnetosphere  and  ionosphere  at  night  to  prevent  F- 
region  decay. 


•  TIGCM  parameterized  inputs  appear  to  be  too  large  at 
times  and  especially  in  later  phase  of  storm  period. 


•  Need  to  continue  comparison  studies  for  different 
periods  using  DE  data,  AE  data  and  data  from  CEDAR 
campaigns. 
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TABLE  2 


Solir  Minimum  (76200)  Sol&r  Miximum  (70050) 


No. 

AT  (5..) 

x  (s^) 

AT  (5*,) 

l 

303.78 

53.5 

303.78 

92.0 

2 

850-900 

20.3 

850-900 

51.7 

3 

150-200 

19.4 

150-200 

47.6 

4 

584.33 

15.1 

300-350 

43.9 

5 

900-950 

14.2 

584.33 

41.7 

6 

629.73 

12.2 

900-950 

34.5 

7 

200-250 

9.9 

200-250 

33.6 

8 

250-300 
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THE  TIGCM  ON  THE  NCAR  CRAY  XMP-48 


PROGRAM  REQUIREMENTS: 

PROGRAM  LENGTH  IN  CORE: 

DISC  REQUIREMENTS: 

Starting  history  volume  (released  after  reading) 
Input  dataset  (released  after  reading) 

History  volume  (written  during  course  of  model  run) 

Copy  of  history  volume  (for  disposal  to  mass  store 
and  save  on  disks) 


.78  megawords 

<  20  mega  words 
0.48  megawords 

<  20  megawords 
(up  to  13  histories) 

<  20  megawords 


TOTAL  DISK  REQUIREMENTS: 


60.48  megawords 


SSD  (OR  DISKS)  REQUIREMENTS: 

Current  and  Updated  Copies  of  Model  Fields 

2  x  1576317  =  3.2  megawords 

Copy  of  Input  Data  Set  .41  megawords 

These  data  sets  are  used  during  time  stepping  and 
employ  optimized  unblocked  data  flow. 

TOTAL  SSD  NEEDS:  3.6  megawords 


•  Model  requires  1722  seconds  of  CPU  time  to  simulate  one  day  using  a  5 
minute  time  step  and  writing  histories  every  hour. 

•  If  a  4  minute  time  step  is  needed,  the  time  for  a  day’s  simulation  increases 
to  2152  secs. 

•  The  model  is  vectorized  and  runs  at  117  megaflops/second  when  time  step¬ 
ping. 

•  A  value  of  50  megaflops/second  is  considered  to  be  reasonably  efficient. 

•  The  code  does  not  make  use  of  multiprocessing. 
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Exreutive  Summary. 

The  Air  Force  has  requirements  for  accurate  knowledge  and  specification  of  the 
state  variables  of  the  Earth's  thermosphere.  Work  done  under  this  contract  will  lead  to  a 
new  model  of  global  thermospheric  density  that  can  be  used  to  specify  and  forecast  neutral 
densities,  temperatures,  and  winds  for  a  wide  range  of  solar  and  geomagnetic  activity.  The 
model  will  be  based  on  simulations  made  with  the  National  Center  for  Atmospheric 
Research  (NCAR)  Thermospheric/Ionospheric  General  Circulation  Model  (TIGCM)  and  on 
data.  The  model  will  be  capable  of  using  real-time  geophysical  indices  or  data  from  ground- 
based  and  satellite  inputs  and  will  provide  neutral  variables  at  specified  locations  and  times 
in  the  altitude  range  90  -  1500  km.  The  model  will  be  semi-empirical  and  will  be  based  on  a 
new  Vector  Spherical  Harmonic  (VSH)  analysis  technique  developed  at  fhe  University  of 
Michigan  that  permits  the  incorporation  of  the  TIGCM  outputs  into  the  model.  The  VSH 
model  will  be  provided  for  use  by  the  Air  Force  in  an  operational  setting  as  a  more  accurate 
version  of  existing  models  of  the  neutral  upper  atmosphere. 

During  the  past  reporting  period,  the  effort  at  the  University  of  Michigan  has 
focused  on  the  further  development  and  improvement  of  the  VSH  model  and  the  design  and 
implementation  of  an  objective  analysis  scheme  for  the  incorporation  of  density 
measurements.  In  addition,  a  research  and  development  science  team  has  been  established 
at  the  University  of  Michigan  to  aid  Prof.  Killeen  in  the  administration  and  implementation 
of  the  model  development  effort.  This  team,  comprising  Prof.  Killeen,  Drs.  Johnson,  and 
Bums  and  Mr.  Raskin,  meets  on  a  regular  weekly  basis. 

Prior  to  the  last  reporting  period,  a  preliminary  version  of  the  VSH  model  was  sent 
to  GL  for  evaluation,  complete  with  12  sets  of  coefficients,  documentation  and  source 
code.  An  initial  response  from  DOAS  (Capt  Osedacz)  was  subsequently  received  (dated  10 
October  1989)  with  helpful  comments.  These  comments  from  the  end  user  of  the  model 
have  been  considered  and  the  suggested  work  related  to  additional  documentation  and  code 
modifications  has  beer,  carried  out  during  the  past  reporting  period.  A  preliminary  version 
of  the  objective  analysis  scheme  for  the  ingestion  of  density  data  has  been  developed  and 
tested,  using  density  measurements  from  the  SETA  instrument  as  a  first  test  case.  Work 
has  also  been  carried  out  to  refine  the  operational  aspects  of  the  model  (i.e.  structured  code 
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development,  documentation,  flags  and  inputs,  data  ingestion  module,  coefficient 
modification  module,  etc.). 

Detailed  sensitivity  and  model  validation  studies  using  the  TTGCM  have  also  been 
carried  out  during  the  past  reporting  period.  The  relevant  sensitivities  of  the  recovered 
densities  to  variations  in  the  different  model  inputs  have  been  quantified.  The  T1GCM  for 
the  Phase  1  Model  has  been  frozen  and  a  significant  number  (18)  model  runs, 
corresponding  to  different  geophysical  conditions,  have  been  performed  and  incorporated 
into  the  VSH  library. 

Scientists  from  the  University  of  Michigan  and  from  NCAR  attended  the  second 
quarterly  review  meeting  in  Boulder  on  21  September  1989.  This  meeting  was  also 
attended  by  many  other  scientists  from  various  institutions  as  well  as  Air  Force  Personnel. 

Plans  for  the  upcoming  reporting  period  include  the  further  development  of  the 
VSH  approach,  leading  to  delivery  of  the  Phase  1  Model  at  the  end  of  FY  1990.  Additional 
key  validation  and  model  improvement  efforts  will  be  carried  ?ut  at  NCAR.  Further 
TIGCM  runs  will  be  made  and  continuing  model  validation  tests  on  these  runs  will  be 
carried  out.  An  improved  treatment  for  the  description  of  the  output  fields  in  the  vertical 
dimension  will  be  designed  and  implemented.  A  theoretical  model  for  neutral  exospheric 
hydrogen  will  be  chosen  and  initial  efforts  undertaken  to  incorporate  it  into  the  VSH.  A 
scheme  for  the  use  of  the  MSIS  model  within  the  VSH  framework  will  be  designed.  MSIS 
predictions  are  useful  to  extend  the  TIGCM  calculation  beyond  the  model  boundaries  and  to 
assist  in  normalization  for  total  densities  in  the  absence  of  experimental  data.  A  second, 
updated  (but  still  preliminary)  version  of  the  VSH  model  will  be  delivered  to  GL  for  further 
evaluation.  Lastly,  a  significant  effort  to  coordinate  the  activities  of  the  neutral  and 
ionospheric  modeling  groups  will  be  made  during  the  upcoming  period. 
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3.1  Progress  during  the  current  reporting  period. 


Progress  prior  to  the  current  reporting  period  was  presented  at  the  second  quarterly 
meeting  in  Boulder,  Colorado  on  September  21  1989.  Progress  during  the  current 
reporting  period  was  be  presented  at  the  third  quarterly  review  meeting  in  Colorado  Springs 
on  January  30,  1990.  The  detailed  agenda  for  this  meeting  is  shown  in  Appendix  3.1. 
Presentations  were  made  by  Prof.  Killeen,  Drs.  Johnson  and  Bums,  and  Mr.  Raskii  of  the 
Unviersity  of  Michigan  and  Dr.  Roble  of  NCAR.  Materials  from  their  presentations  are 
included  as  appendicies  to  this  report 

Project  personnel  and  their  responsibilities  are  shown  in  Appendix  3.1.  Appendix  3.2 
reviews  project  milestones  through  the  5-year  contract  and  summarizes  the  current  status  of 
the  project.  In  spite  of  a  late  start  for  the  first  fiscal  year  of  the  contract  in  March,  1989, 
progress  in  model  development  is  on  schedule.  Tasks  1  through  3  have  been  completed 
and  the  science  team  is  currently  continuing  work  on  tasks  4  through  7  to  achieve  an  on 
time  delievery  of  the  Phase  1  Model  at  the  end  of  the  second  fiscal  year. 

The  status  of  the  model  and  plans  for  its  further  development  and  refinement  are  presented 
in  Appendix  3.3.  A  plan  for  data  flow  and  manipulation  within  the  VSH  model  has  been 
established,  including  the  capability  for  modification  of  both  coefficient  files  and  VSH 
output  fields.  Default  density  models  (MSIS,  exospheric  H)  will  be  utilized  at  altitudes 
outside  the  range  of  TIGCM  validity. 

Total  neutral  density,  which  is  calculated  from  the  mass  mixing  ratios,  has  been  added  as 
an  output  field.  We  are  currently  considering  adding  an  additional  set  of  coefficients  for  the 
total  density  to  the  coefficient  files  to  expedite  potential  objective  analysis  schemes. 

The  mathematical  formulation  for  the  several  transformations  utilized  in  the  VSH  model  are 
specified  in  Appendix  3.3.  Alternatives  to  the  Fourier  time  representation,  such  as 
Chebyshev  polynomials,  are  currently  being  considered  for  application  to  temporally 
variable  storm-time  conditions.  In  order  to  adequately  describe  high  latitude  and  low 
altitude  fields,  we  are  also  investigating  alternatives  to  the  spatial  transform  currently 
implemented.  As  more  altitudes  are  retained  from  the  TIGCM,  mathematical 
representations  for  the  vertical  variation  of  the  coefficient  fields  must  be  developed.  We  are 
currently  investigating  potential  formulations. 
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An  output  field  modifier  has  been  developed  for  the  mass  density.  We  also  plan  to  develop 
an  operational  coefficient  modification  module  (currently  coefficient  files  can  only  be 
modified  on  the  NCAR  Cray). 

Our  responses  to  the  10  October  1989  memo  from  Capt  Osedacz  are  indicated  in  Appendix 
3.3.  Specifically,  we  have  completed  the  addition  of  mass  density  as  an  output  field.  We 
have  provided  coefficient  files  in  ASCII  format,  and  additionally  provided  a  conversion 
program  back  to  binary.  All  library  subroutines  have  been  documented  internally,  and  the 
units  of  all  output  variables  have  been  specified.  We  are  in  the  process  of  developing  a 
technical  manual  for  the  VSH  model. 

The  status  of  NCAR  TIGCM  runs  and  coefficient  file  production  for  Phase  1  Model 
development  is  summarized  in  Appendix  3.4.  Presently  18  TIGCM  Phase  1  runs  have 
been  completed  for  solar  minimum  and  maximum  conditions,  three  levels  of  geomagnetic 
activity,  and  solstices  and  equinox.  To  date,  eleven  of  the  corresponding  coefficient  files 
have  been  generated. 

The  objective  analysis  schemes  utilized  in  the  VSH  mode)  are  reviewed  in  Appendix  3.5. 
The  objective  analysis  scheme  for  neutral  wind  data  (Dynamics  Explorer  measurements) 
was  developed  in  the  last  reporting  period.  This  period,  we  have  developed  a  preliminary 
objective  analysis  technique  for  ingesting  total  density  measurements.  In  this  scheme,  only 
the  variation  of  density  with  altitude  is  considered,  and  the  VSH  output  fields  are  adjusted 
to  match  the  measured  altitude  profile.  We  are  currently  using  linear  fits  to  the  measured 
and  model  profiles  in  this  transformation.  Two  examples  of  the  capability  of  this  technique 
are  shown  in  Appendix  3.5. 

Potential  improvements  to  this  preliminary  scheme  include  the  utilization  of  more 
appropriate  functions,  such  as  quadratic  or  higher  order  forms.  These  will  certainly 
become  necessary  as  the  altitude  range  of  the  available  data  increases.  Additionally,  an 
option  to  modify  coefficient  files  directly  will  provide  computational  efficiency.  A  default 
to  TIGCM  adjusted  to  match  MSIS  absolute  density  at  low  altitude  and  to  MSIS  augmented 
by  exospheric  H  above  the  maximum  TIGCM  altitude  will  also  be  necessary.  Techniques 
for  density  field  modification  including  latitudinal,  longitudinal,  and  universal  time 
dependences  need  to  be  developed.  Finally,  the  region  of  influence  for  perturbations  on 
average  profiles  will  need  to  be  considered. 
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TIGCM  developments  and  the  results  of  the  Task  3  sensitivity  study  are  summarized  in 
Appendix  3.6.  There  has  been  considerable  progress  in  the  specification  of  physical 
processes  in  the  TIGCM.  Among  other  improvements,  viscous  heating  and  0+  plasma 
exchange  have  been  added,  the  tidal  specification  has  been  improved,  and  altitude 
dependent  gravity  has  been  added.  Sensitivity  studies  for  several  variables  have  been 
performed  which  show  the  relative  change  in  TIGCM  temperature  and  density  to  10% 
changes  in  solar  EUV  and  SRC,  ±  10  nT  changes  in  By,  ±  10  kV  changes  in  potential 
drop,  and  changes  in  the  lower  boundary  tidal  specification. 


3 2  Plans  for  the  next  reporting  period. 

Plans  for  the  next  quarter  of  effort  relate  to  the  work  statement  given  in  Appendix 
3.2.  Specifically,  work  on  tasks  4-7  will  continue.  Major  action  items  for  the  next  3 
months  are  indicated  in  Appendix  3.7.  These  include  incorporation  of  an  exospheric 
hydrogen  model  above  500  km,  an  MSIS  default  scheme,  and  continued  development  and 
implementation  of  the  library  scheme  for  Phase  1  Model  delivery.  We  will  also  deliever  a 
second  preliminary  test  version  of  the  VSH  to  the  Space  Forecast  Center  for  continuing 
evaluation,  and  develop  further  coordination  of  efforts  with  the  ionospheric  modelling 
group.  Continuing  work  on  Tasks  4  through  7  will  include  further  development  and 
refinement  of  the  VSH  algorithm,  specifically  the  vertial  and  storm  time  representations  to 
be  utilized.  Refinement  and  development  of  field  and  coefficient  modifiers  will  continue. 
The  TIGCM  sensitivity  study  will  also  continue,  as  will  Phase  1  TIGCM  runs  for  model 
delievery  at  the  end  of  FY  1990. 
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Appendix  3.1 
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VSH  Model  Development 
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VSH  -  AFGUAWS,  Colorado  Springs  1/28/90  RJ  9 
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Review  of  Schedule 
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VSH  -  AFGL/AWS,  Colorado  Springs  1/28/90 
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VSH  -  AFGL/AWS,  Colorado  Springs  1/28/90  rj  4 
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VSH  -  AFGUAWS,  Colorado  Springs  1/28/90 


VSH  Model  Development 


VSH  -  AFGL/AWS,  Colorado  Springs  1/28/90 
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Mich/NCAR 
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VSH  -  AFGUAWS,  Colorado  Springs  1/28/90  BJ  7 
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VSH  -  AFGL/AWS,  Colorado  Springs  1/28/90  RJ  10 
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VSH  -  AFGUAWS,  Colorado  Springs  1/28/90 


VSH  -  AFGL/AWS,  Colorado  Springs  1/28/90 
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VSH  -  AFGL/AWS,  Colorado  Springs  1/29/90 
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VSH  DATA  FLOW  DIAGRAM 


VSH-AFGL/AWS,  Colorado  Springs  V2M0 
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VSH-AFGL/AWS,  Colorado 
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VSH-AFGUAWS,  Colorado  Springs  1/29/90 
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VSH-AFGUAWS,  CoWx*$o  Spring*  tWSO 
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VSH-AFGL/AWS,  Colorado  Springs  1/29/90 
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POLYNOMIAL  INTERPOLATION  OVER  ALTITUDE 


TIME  TRANSFORM 
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VSH-AFGL/AWS,  Colorado  Spring*  I/S' 
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VSH-AFGUAWS,  Colorado  Spring*  VZW90 
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•  Truncation  level  is  N=25  (meridional  resolution), 

M=  5  (zonal  resolution). 

•  The  750  coefficients  (per  variable)  are  reduced  to  3. 


VECTOR  SPATIAL  TRANSFORM 
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VSH-AFGLVAWS,  Cotoraco  Spring*  1/29 /VO 
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jH-AFGUAWS,  Colorado  Springs  1/29/90  Mil 


NGOING  DEVELOPMENT  EFFORTS 
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-  Empirical  Orthogonal  Functions  (EOPs) 


EXAMPLE  OF  TIDAL  WINDS 
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VSH-AFGL/AWS,  Colorado  Springs  1/29/90 


VSH-AFGUAWS,  Colorado  Springs  1/29/90 
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VSH-AFGUAWS,  Colorado  Springs  1/29/90  HR16 
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VSH-AFGUAW3,  Colorado  Spring* 


VSH  MODEL  ACTION  ITEMS 
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VSH-AFGL/AWS,  Col"  ado  Springs  1/29/90  RR18 


103 


Write  technical  manual  of  VSH  model 


104 


Run  Status 
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VSH  -  AFGL/AWS,  Colorado  Springs  1/28/90 
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VSH  -  AFGL/AWS,  Colorado  Springs  1/28/90  AB  2 


Selection  Criteria 
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VSH  -  AFGL/AWS,  Colorado  Springs  1/28/90  AB  3 


General 


VSH  •  AFGUAWS,  Colorado  Springs  1/28/90 
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VSH  •  AFGUAWS,  Colorado  Springs  1/28/90 
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VSH  -  AFGL/AWS,  Colorado  Springs  1/28/90 


OBJECTIVE  ANALYSIS  SCHEME 
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Use  regular  vector  spherical  harmonic  truncation  scheme 


DE-2  NEUTRAL  WIND 
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•  3 .0  <  Kp  <  6.0,  By  <  0.0  #  of  Orbits  =  1 26 


Comparison  of  Averaged  (1983-88)  Sondrestrom 
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RJ  9 


VSH  -  AFGUAWS,  Colorado  Springs  1/28/90 


Objective  Analysis  of 


VSH  -  AFGL/AWS,  Colorado  Springs  1/28/90 
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VSH  -  AFGUAWS.  Colorado  Springs  1/28Z90  RJ  18 


Objective  Analysis  of 
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VSH  -  AFGL/AWS,  Colorado  Springs  1/20/90  RJ  20 
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TIGCM  Developments 


1.  Solar  EUV/UV  Flux  Model 

•  Based  on  F107  and  F10  7  (Hinteregger  data) 

•  Based  on  Lyman  -a  and  Fj07  (SME  data  and  Hinteregger  data) 

•  Compared  with  recent  rocket  measurements 

•  Largest  uncertainties  at  X  <  20  nm  (important  for  E-region,  ne 
and  NO) 

2.  Improved  Auroral  Model 

•  Evans  and  Fuller-Rowell  particle  precipitation 

•  Convection  model 

a)  Heelis,  analytical 

b)  Maynard  and  Rich,  analytical 

c)  Emery  modifications  for  time-dependent  forcings 

d)  AMIE 

•  Needs  considerable  development  for  time-dependent  forcings 

3.  Magnetosphere/Ionosphere  Plasma  and  Energy  Exchange 

•  Considerable  improvements  achieved  by  specifying  fluxes 

•  Empirical  parameterizations  used 

•  Need  to  develop  physical/interactive  model 

4.  Tidal  Boundary  (Semi-Diurnal  Tides) 

•  Amplitude  and  phases  from  Vial/Forbes  model  for  each  month 
of  the  year 

•  Verified  with  data  (except  September  tides) 
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TIGCM  Developments  (Cont.) 


5.  Introduced  1-1  Tide  into  TIGCM 

•  Important  for  lower  thermosphere 

•  Amplitude  and  phase  not  well  defined 

6.  Included  Viscous  Heating  into  TIGCM  Thermodynamic  Equation 

1-2%  Changes  in  Density) 

7.  Improve  Specification  of  Height  of  Constant  Pressure  Surface 

•  Introduce  variable  gravity  (~  30  km  at  upper  B.C.) 

•  TIGCM  densities  were  previously  high  relative  to  data  and 
MSIS-86 

•  Much  better  agreement  with  variable  gravity 

•  Future  improvements  require  connection  with  lower  atmosphere 
(model  atmosphere  to  specify  heights) 

•  Normalization  is  improved  but  needs  further  development 

8.  Improve  TIGCM  Processors 

•  Comparisons  with  data  and  MSIS-86 

•  Analysis  of  physical  processes 

•  Rapid  visualizations 
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TIGCM  Validation 


A  great  deal  of  progress  has  been  made  in  substantiating  the  TIGCM  per¬ 
formance  agamst  real  data  and  also  in  documenting  TIGCM  performance 
versus  MSIS-86  model  predictions.  At  least  five  papers  are  in  various  stages 
of  publication  documenting  TIGCM  performance.  In  addition,  the  model  is 
being  used  by  at  least  eight  d  liferent  scientific  collaborators  to  verify  model 
p.  rformance  by  comparison  with  various  ground  based  measurements 
which  include: 

a)  Effect  of  tides  from  the  lower  atmosphere  on  thermospheric  den¬ 
sities  --  Fesen,  U.  o'  Colorado;  Forbes,  Boston  V. 

b)  Ionospheric  storm  perturbations  -  Codrescu,  Boston  U.;  Crow¬ 
ley,  Lowell  U. 

c)  Airglow  emissions  -  Solomon,  NCAR 

d)  Thermospheric  density  and  composition  --  Burns  and  Killeen,  U. 
of  Michigan 

e)  Vorticity  and  divergent  flows  -  Thayer  and  Killeen,  U.  of 
Michigan 

f)  Nitric  oxide  --  Fesen,  U.  of  Colorado;  Gerard.  U.  of  Leige 

g)  South  polar  cap  dynamics  --  Smith,  U.  of  Alaska;  Hernandez.  U. 
of  Washington 

h)  Radar  derived  winds  -  Wickwar,  Utah  State;  Hagan,  Millstone 
Hill;  Burnside,  Arecibo  Obs. 

All  of  these  studies  are  providing  inputs  to  help  validate  model  performance 
and  suggest  areas  where  further  TIGCM  improvements  are  necessary  to 
achieve  the  Air  Force  overall  density  goak.  The  specific  case  studies  that 
are  underway  include: 

a)  March  21  -  April  10,  1979  —  SETA,  radar  data 

b)  January  21-31,  1983  -  Comprehensive  DE  measurements 

c)  February  -4-9,  1986  major  storrn  -  Aerospace  airglow  and  radar 
measurements 

d)  September  20-28,  198-4  --  ETS  study 

e)  March  13-14,  1989  major  storm 

f)  LTCS,  September  21-25,  1987,  December  5-10,  1988  period 
plus  there  are  others  associated  with  CEDAR,  WITS,  and  DE. 


TIGCM  Studies  in  Progress 


1.  TIGCM  Sensitivity  Studies 

a)  Increase  solar  EUV  by  10% 

b)  Decrease  solar  EUV  by  10% 

c)  Change  By  from  0  to  +10 

d)  Change  By  from  0  to  -10 

e)  Change  <f>  from  40  to  50  kV 

f)  Change  from  40  to  30  kV 

1)  Solar  min;  March  equinox,  December  solstice 

F 10  7  =  67,  F10  7  =  72,  P  =  5  GW,  4>  =  40  kV 

2)  Solar  max;  December  solstice 

F107  =  243,  F107  =  215,  P  =  11  GW,  $  =  60  kV 

g)  Increase  solar  SRC  by  10% 

h)  Change  tidal  specification 

•  Results  are  presented  as  %  difference  from  unperturbed  run 

2.  TIGCM/MSIS-86  Comparisons 

•  18  TIGCM  cases  have  been  run  for  density  model  development 

•  Develop  normalization  scheme 

3.  Continue  March  1979  Study  using  SETA  Data 

•  20  day  run,  March  20  -  April  10,  1979 

•  Percentage  errors  from  SETA  measurements 

a)  TIGCM 

b)  MSIS-86 

4.  DE  Satellite  Study,  January  20-30,  1983 

•  Tn  ,  Vn  ,  nn  ,  nt  ,  Te  ,  comparisons  along  satellite  track 
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I  HAVE  MADE  TIGCM  RUNS  FOR  A  CONSISTENT  SET  TO  BE  USED  IN  AIR  FORCE 
DENSITY  SPECIFICATION  MODEL  DEVELOPMENT 

A.  SOLSTICE  SOLAR  MINIMUM  OA Y « 76355 , F 107 =  67 . F 1 07 A = 72 

OATA  *ROBL  E / RGR88 /D VWSOL 

1 .  P*5GW.CTPOTEN=30KV 

HI  STORY *R0BLE/RGR90/SQLN0l ,02  FROM  1.0,0  TO  2,0,0  HOURLY 

2 .  P* 1 1 GW . CTP0TEN*6OKV 

HI  STORY =R0BL£/RGR90/SQLN1 1 . 12  FROM  1,0,0  TO  2.0,0  HOURLY 

3 .  P=33GW . CTP0TEN=9OKV 

HI  STORY *ROBL£ /RGR90/SOLN2 1 ,22  FROM  1.0,0  TO  2,0.0  HOURLY 

B.  SOLSTICE  SOLAR  MAXIMUM  DAY =79355 . F 107 =243 . F 107A=2 15 

1 . P=5GW,CTP0TEN=3OKV 

HISTORY=ROBLE/RGR90/SOLX01 .02  FROM  1.0,0  TO  2,0,0  HOURLY 

2 .  P= 1 1GW ,CTP0TEN=6OKV 

HISTORY=ROBLE/RGR90/SOLX 1 1 , 12  FROM  1,0,0  TO  2,0.0  HOURLY 

3 .  P=33GW , CTP0TEN=90KV 

HI  STORY = R08LE /RGR9O/S0LX2 1 ,22  FROM  1.0,0  TO  2.0.C  HOURLY 
C  EQUINOX  SLOAR  MINIMUM  DA Y = 76080 . F 107 =67 . F 107A = 72 . 

1 .  P  =  5GW . CTP0TEN=30KV 

HI  STORY *R0BLE/RGR90/EQSN0 1 ,02  FROM  1,0,0  TO  2,0,0  HOURLY 

2. P* 1 1GW.CTP0TEN=60KV 

HI  STORY =ROBLE/RGR9O/E0SN1 1 , 12  FROM  1.0,0  TO  2.0.0  HOURLY 

3 .  P  =  33GW , CTPOTEN*90KV 

HIST0RY=R0BLE/RGR90/E0SN21 ,22  FROM  1,0.0  TO  2,0.0  HOURLY 

D.  EOUINOX  SOLAR  MAXIMUM  DAY =79080 , F 107 = 243 . F 107 A =2 15 

1 .P«5GW , CTP0TEN=3OKV 

HISTORY  =ROBLE /RGR90/ EQSX03 , 04  FROM  1,0,0  TO  2,0,0  HOURLY 
2  - P* 1 1GW.CTPOTEN=60KV 

HI  STORY  =  ROBLE /RGR90/ EQSX 1 1 ,  12  FROM  2.0,0  TO  3.0.0  HOURLY 
3 . P= 1 1GW.CTPOTEN=90KV 

HISTORY  =RQBL E /RGR90/ EQSX2 1 , 22  FROM  1,0.0  TO  2,0,0  HOURLY 

E.  SUMMER  SOLSTICE  SOLAR  MINIMUM  DA Y « 76 1 82 , F 107 =67 , F 107A =72 

1 .  P=5GW.CTP0TEN=3OKV 

HISTORY  =  R0BLE/RGR90/SUMN05 , 06  FROM  2.0,0  TO  3,0,0  HOURLY 

2 .  P= 1 1GW.CTP0T  EN  =  60KV 

HI  STOR Y = ROB L E /RGR90/SUMN 1 3 ,  14  FROM  3,0.0  TO  4.0.0  HOURLY 

3 .  P*  33GW . CTP0TEN  =  9OKV 

HIST0RY=R0BLE / RGR90/ SUMN2 3,24  FROM  3.0,0  TO  4,0,0  HOURLY 

F.  SUMMER  SOLSTICE  SOLAR  MAXIMUM  D A Y = 7 9 1 8 2 , F 1 07 = 2 4 3 , F 1 07 A = 2 1 5 

1 .  P  =  5GW , CTPOTEN  =  30KV 

HI STOR Y = ROBL E /RGR90/ SUMXO 1 ,02  FROM  1.0.0  TO  2,0.0  HOURLY 

2 .  P= 1 1  GW , CTP0TEN=6OKV 

H I STOR Y = ROB L E /RGR90/ SUMX 1 1 . 12  FROM  1,0,0  TO  2,0.0  HOURLY 

3 .  P=33GW . CTP0TEN=90KV 

HIST0RY=R08LE /RGR90/SUMX2 1.22  FROM  1,0.0  TO  2.0.0  HOURLY 
(NOTE  THAT  THE  NEW  MODE L/RGR90/H I STOR I E S  ARE  UT  =MODE  L  T I  ME ( F I NAL L Y ) 
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Critical  Inputs  Required  during  Quiet  Geomagnetic  Activity 


Thermosphere/Ionosphere  System 


Zero  Order  First  Order 


Higher  Order 


Solar  EUV/UV 
5  <  X  <  200  nm 

Particle  Precipitation 
(GW) 

Cross-Cap  Potential 
(kV) 

Bz,  By,  Configuration 
of  Potential  Pattern 

Tidal  Amplitudes 
and  Phases 

Density  Measurements 
(normalization) 

Magnetosphere/Ionosphere 
Plasma  and  Energy 
Exchange 
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Critical  Inputs  Required  during  Disturbed  Geomagnetic  Activity 


Thermosphere/Ionosphere  System 


Zero  Order _ First  Order _ Higher  Order 

Solar  EUV/ITV 
5  <  X  <  200  nm 


Bz,  By,  Configuration 
of  Potential  Pattern 

Magnetosphere/Ionosphere 
Plasma  and  Energy 
Exchange 

Tidal  Amplitudes 
and  Phases 

Density  Measurements 
(normalization) 


Particle  Precipitation 
(GW) 

Cross-Cap  Potential 
(kV) 
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MINIMUM  598  2,  MAXIMUM  838  6,  CONTOUR  INTERVAL  15 


MINIMUM  0  9060*10  'V  MAXIMUM  1  6000*10  '*.  CONTOUR  INTERVAL 


MINIMUM  6  789,  MAXIMUM  10  505,  CONTOUR  INTERVAL 


MINIMUM  1  017.  MAXIMUM  2  740,  CONTOUR  INTERVAL 


MINIMUM  -10  193,  MAXIMUM  -6.588.  CONTOUR  INTERVAL 
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MINIMUM  -9.00,  MAXIMUM  9.81,  CONTOUR  INTERVAL  1 
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MINIMUM  282.  MAXIMUM  3.959.  CONTOUR  INTERVAL  20 


T1GCM  NEUTRAL  TEMPERATURE  (DEG  K)  (Difference) 
%D1FF  (Potential  =  50kV)  UT=12.0  HT  =  300.0 


MINIMUM  MAXIMUM  (J ‘Jby  CONTOUR  INTERVAL 


DENSITY 


MINIMUM  -2  397,  MAXIMUM  7  979,  CONTOUR  INTERVAL 
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Summary  and  Conclusions 


•  Considerable  Improvement  in  TIGCM  Physical  Processes 

a)  Viscous  heating 

b)  0+  plasma  exchange 

c)  Tidal  specifications 

d)  etc. 

•  Performed  a  Large  Number  of  Sensitivity  Studies 

a)  Solar  EUV  variations 

b)  Solar  SRC  variations 

c)  By  and  4>  variations 

d)  Tides 

e)  Others 

•  Improved  TIGCM  Normalization;  Better  Agreement  with  Data 

•  Continued  with  Data  Validation  Studies;  Significant  TIGCM 
Improvements 
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Major  Action  Items  For 
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S,  Colorado  Springs  1/29/90 
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VSH  -  AFGL/AWS,  Colorado  Springs  1/29/90  RJ  22 
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Executive  Summary 

The  Air  Force  has  requirements  for  accurate  knowledge  and  specification  of  the 
state  variables  of  the  Earth's  thermosphere.  Work  done  under  this  contract  will  lead  to  a 
new  model  of  global  thermospheric  density  that  can  be  used  to  specify  and  forecast  neutral 
densities,  temperatures,  and  winds  for  a  wide  range  of  solar  and  geomagnetic  activity.  The 
model  will  be  based  on  simulations  made  with  the  National  Center  for  Atmospheric 
Research  (NCAR)  Thermosphere/Ionosphere  General  Circulation  Model  (T1GCM)  and  on 
data.  It  will  be  capable  of  using  real-time  geophysical  indices  or  data  from  ground-based 
and  satellite  inputs  and  will  provide  neutral  variables  at  specified  locations  and  times  in  the 
altitude  range  90  -  1500  km.  This  hybrid  theoretical  and  semi-empirical  model  will  be  based 
on  a  new  Vector  Spherical  Harmonic  (VSH)  analysis  technique  developed  at  the  University' 
of  Michigan  that  permits  the  incorporation  of  the  TIGCM  outputs  and  data  into  the  model. 
The  VSH  model  will  be  provided  for  use  by  the  Air  Force  in  an  operational  setting  and  will 
be  a  more  accurate  version  of  existing  models  of  the  neutral  upper  atmosphere. 

During  the  reporting  period,  the  effort  at  the  University  of  Michigan  has  focused  on 
the  further  development  and  improvement  of  the  VSH  model.  The  MSIS  empirical  model 
has  been  incorporated  into  the  VSH,  extending  the  range  of  the  model  above  the  upper  limit 
of  the  TIGCM  to  1500  km.  A  scheme  has  been  developed  by  which  the  model  interpolates 
between  values  in  the  spectral  coefficient  library  to  user-specified  geophysical  conditions. 
Additionally,  several  changes  in  VSH  variables  and  outputs  have  been  made.  Work  on 
extending  the  TIGCM  run  set  for  Phase  1  Model  delivery  has  continued.  Eighteen 
coefficient  files  are  pi^ntlj  useo  io  run  me  eating  VSH  model,  incorporating  changes  in 
season,  changes  in  solar  cycle  and  changes  in  geomagnetic  activity.  Additionally,  several 
changes  in  VSH  variables  and  outputs  have  been  made. 

Plans  for  the  upcoming  reporting  period  include  continued  development  of  the  VSH 
approach,  leading  to  delivery  of  the  Phase  1  Model  at  the  end  of  FY  1990.  Additional  key 
validation  and  model  improvement  efforts  will  be  carried  out  at  NCAR.  Further  TIGCM 
runs  will  be  made  and  continuing  model  validation  tests  on  these  runs  will  be  carried  out. 
An  improved  treatment  for  the  description  of  the  output  fields  in  the  vertical  dimension  will 
be  designed  and  implemented.  Work  on  data  ingestion  techniques  will  continue.  An 
improved  theoretical  model  for  neutral  exospheric  hydrogen  will  be  chosen  and  initial 
efforts  undertaken  to  incorporate  it  into  the  VSH.  Lastly,  an  effort  to  coordinate  the 
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activities  of  the  neutral  and  ionospheric  modeling  groups  will  continue  during  the  upcoming 
period. 


4.1  Progress  during  the  current  reporting  period. 

^  Appendix  4.1  reviews  project  milestones  through  the  5-year  contract  and 

summarizes  the  current  status  of  the  project.  In  spite  of  a  late  start  for  the  first  fiscal  year  of 
the  contract  in  March,  1989,  progress  in  model  development  is  on  schedule.  Tasks  1 
through  3  have  been  completed  and  the  science  team  is  currently  continuing  work  on  tasks 
4  through  7  to  achieve  an  on  time  delivery  of  the  Phase  1  Model  at  the  end  of  the  second 
fiscal  year. 

The  MSIS  empirical  model  has  been  incorporated  within  the  VSH,  extending  the 
nominal  range  of  the  model  to  1500  km.  For  each  of  the  TIGCM  runs  used  as  the  base  set 
for  the  VSH,  the  mid-latitude  average  mass  density  is  compared  to  the  corresponding 
average  from  MSIS.  Their  relative  magnitudes  define  a  normalization  factor  that  is  applied 
as  an  output  field  modifier  for  mass  density  at  the  last  step  in  the  density  calculation.  As 
part  of  the  MSIS  model  implementation,  hydrogen  density  has  been  added  as  a  13th  output 
variable.  At  any  altitude,  hydrogen  densities  are  obtained  from  the  MSIS  model.  Above 
the  upper  altitude  limit  of  the  TIGCM,  MSIS  is  called  to  specify  all  constituent  densities. 

Work  has  progressed  on  the  design  of  the  coefficient  file  library  and  its  utilization  in 
the  VSH  model.  Rather  than  specifying  a  TIGCM  run,  the  user  specifies  F10.7,  Julian 
day,  and  Kp  values.  The  VSH  model  then  uses  these  inputs  to  determine  appropriate 
spectral  coefficients  through  linear  combination  of  those  determined  from  individual 
TIGCM  runs  (in  the  coefficient  file  library).  Kp  is  converted  to  hemispheric  power  for 
purposes  of  finding  the  linear  combination  weights  over  geomagnetic  activity.  The  weights 
for  the  day  of  year  interpolation  are  found  via  a  trigonometric  scheme. 

Additionally,  several  other  changes  have  been  made  to  the  VSH.  Meridional  winds 
are  now  specified  as  positive  northwards,  rather  than  positive  southwards.  Colatitude  has 
been  replaced  by  latitude.  Finally,  mass  density  has  been  added  as  a  14th  output  variable. 

*  We  are  nearing  completion  of  a  technical  manual  for  the  VSH  model.  A  second  preliminary 

version  of  the  VSH,  updated  with  changes  reported  in  the  last  quarterly  report,  was 
delivered  to  GL  for  further  evaluation. 
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Presently  18  TIGCM  Phase  1  runs  have  been  completed  for  solar  minimum  and 
maximum  conditions,  three  levels  of  geomagnetic  activity,  and  solstices  and  equinox.  All 
of  the  corresponding  VSH  spectral  coefficient  files  have  been  generated.  New  coefficient 
files  are  being  planned  at  present,  so  that  a  more  accurate  vertical  representation  of  the  state 
variables  can  be  obtained,  and  so  that  data  can  be  better  integrated  with  the  model.  Further 
work  is  being  undertaken  to  identify  new  TIGCM  runs  needed. 


4.2  Plans  for  the  next  reporting  period. 

Plans  for  the  next  quarter  of  effort  relate  to  the  work  statement  given  in  Appendix 
4.1.  Specifically,  work  on  tasks  4-7  will  continue.  Major  action  items  for  the  next  quarter 
include  continued  refinement  of  both  the  library  scheme  for  Phase  1  Model  delivery  and  the 
VSH  algorithm,  specifically  the  vertical  and  storm  time  representations  to  be  utilized. 
Development  of  field  and  coefficient  modifiers  and  data  ingestion  techniques  will  continue. 
The  TIGCM  sensitivity  study  will  also  continue,  as  will  Phase  1  TIGCM  runs  for  model 
delivery  at  the  end  of  FY  1990.  We  will  continue  to  develop  further  coordination  of  efforts 
with  the  ionospheric  modelling  group. 
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Model  Development  Progress 

-  Milestone  Completed _ >  -  Work  in  Progress  O  -  Work  not  yet  begun 
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